Modelling the effect of charge noise on the exchange interaction between spins 
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We describe how the eflect of charge noise on a pair of spins coupled via the exchange interaction 
can be calculated by modelling charge fluctuations as a random telegraph noise process using prob- 
ability density functions. We develop analytic expressions for the time dependent superoperator of 
a pair of spins as a function of fluctuation amplitude and rate. We show that the theory can be 
extended to include multiple fluctuators, in particular, spectral distributions of fluctuators. These 
superoperators can be included in time dependent analyses of the state of spin systems designed 
for spintronics or quantum information processing to determine the decohering effects of exchange 
fluctuations. 
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Introduction 

The exchange interaction is of increasing importance 
in the study of controllable quantum mechanics using 
solid-state systems. As well as being fundamentally 
important in manybody physics, it is this interaction 
which is often used to mediate spin flips or entangle- 
ment in spintronics and quantum information processing 
(QIP) [H, i, i, i, i, i, For these reasons there has 
been considerable study recently on the ori gin and con- 
trol of the exchange interaction [1, [^, [l^, [Til Il2l. [l3| . For 
applications involving the time varying control of the ex- 
change interaction, such as QIP, the stability in time of 
this interaction is of crucial importance. As the origin 
of the exchange interaction is essentially the overlap of 
electron wave functions, the interaction strength is sen- 
sitive to the local charge environment. Recent work 
has shown that the dependence of the exchange interac- 
tion is approximately linearly dependent on fluctuations 
in the local electric field. 

In this paper we develop a general framework with 
which the effect of these fluctuations can be analyt- 
ically included in time dependent calculations of the 
state of a sp in system. Using the superoperator formal- 
ism [ill, Ea Ej [H, j ^6 derive an expression for 
a pair of spins as a function of time, depending on the 
exchange fluctuation amplitude and rate. The extension 
to multiple fluctuators, in particular, spectral distribu- 
tions of fluctuators is also considered. The formalism can 
be used to investigate the effect of exchange fluctuations 
on spintronics, quantum control schemes and speciflcally 
quantum error correction (QEC). 



I. THE NOISE MODEL 

We begin by studying the exchange coupling Hamil- 
tonian in the presence of a single charge fluctuator with 
the aim of understanding the decohering effects of the 



fluctuator. The Hamiltonian for the process is 

H (t) = .J (t) a-i ■ (72. (1) 

The exchange coupling J (t) varies in time due to a ran- 
dom telegraph noise (RTN) process, rj (t) and we assume 
a net effect of the form 



J (t) = Jo + mi (t) 



(2) 



where ?/ (t) describes the fluctuator. This RTN process 
couples with strength a (ultimately dependent on the 
distance between the coupled spins and the fluctuator) 
to the bare exchange term, Jg. The time evolution of 
the system can then be described by the density matrix 
master equation 



= [H {t),p{t)]. 



(3) 



where p [t] is the density matrix of the system. Addi- 
tional terms can be added to this master equation to 
also model non-unitary evolution, such as decohering pro- 
cesses. 

As a matter of convenience we may re-express the sys- 
tem evolution in superoperator form. In superoperator 
form the density matrix is given a vector representation, 
denoted by p(t), by transforming the matrix into a single 
column, one row at a time plj . A superoperator P (t) 
contains all the evolution of the system (both unitary 
and non-unitary) 



p it) ^ Pit) Pit). 



(4) 



For purely Hamiltonian evolution the superoperator P (t) 
can be written down in terms of H it) and the identity 
operator 



Fit) = -i Hit)(g)I -I(g>H it)^ 
The superoperator simplifies to 

Pit) = -iJit) (TH, 



(5) 



(6) 
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for the Hamiltonian we consider. Here, cth is the Heisen- 
berg interaction in superoperator form. If the Hamilto- 
nian is time independent, then the superoperator P is 
also time independent and the density matrix at some 
time t is 



Pit) 



Pito) , 
^{t)p{to) 



(7) 
(8) 



given an initial state p{to). We show how this time in- 
dependent formalism is relevant to our problem shortly. 

The RTN process, ri{t), is modelled as in Ref. [23 |. 
The noise fluctuates randomly between -1 and 1 with the 
frequency of the fluctuations controlled by the correlation 
time 1/A. Here, A is the typical frequency of jump times, 
where the jump time instants are, 



' 1 



(9) 



and the pj are random numbers such that pj £ (0,1). 
The noise process r] {t) is described as 



E.e(t-t,) 



(10) 



where Q (t) is the Heaviside step function, and rj (t) can 
fluctuate between ±77(0). We choose |?7(0)| = 1 and 
control the coupling strength via a as in Eq. [21 

The density matrix evolution for our system can be 
found by numerically averaging over many such noise his- 
tories r] (t) to obtain the correct system dynamics. For 
an initial state p{to), 



1 ^ 



(11) 



k=l 



where the {Uk} are the evolution operators for trajecto- 
ries r]k (t) . Since the Hamiltonian (Eq. [T]) commutes with 
itself at all times. 



[H{to),H{t)]=0, 
the {Uk} may be expressed as 

C/fe(t,to) = C/-(t-)C/+(t+), 



(12) 



(13) 



where i_ and describe the total time the fluctuator 
exists in the -1 and -1-1 states respectively for a particular 
noise history, and 



i{Jo±a}(T 1 -Cr 2* 



(14) 



Using the result of Eq. [8] for the superoperator form of a 
density matrix governed by a time independent Hamilto- 
nian we re-express Eq. [11] such that 



1 ^ 

AT— >oo iV 



(15) 



The ensemble averaged superoperator, Q (t), is the aver- 
age of all the individual trajectory superoperators q^, (t), 



1 ^ 



(16) 



k=l 



This implies that Q (t) may be constructed by numeri- 
cally averaging over many noise histories. The averaging 
is crucial in obtaining the correct system dynamics, as the 
RTN is a stochastic process and so there are many unique 
noise trajectories. Averaging over these noise trajectories 
results in non-unitary evolution despite the Hamiltonian 
being strictly unitary. 

Conversely, it is possible to derive Q (t) analytically 
by describing the stochastic RTN using an appropriate 
probability density function (PDF). By considering all 
unique (t) as a function of the average fluctuator state 
£^ = \r]{0)\{t+ -t^) /T, weighted by a PDF giving the 
occurrence likelihood of the average fluctuator state, and 
integrating this over all possible ^, the resulting expres- 
sion for Q (t) is 



Q(0 = / qi{t)n{^,T)d^. 



(17) 



Here, (t) is the unique individual superoperator cor- 
responding to a particular value of ^ and ft (^, T) is the 
PDF, which determines the probability that during the 
time interval T, the average fluctuator state is ^. In sec- 
tion [TT] we show how to specify the PDF, so that we can 
use it to analytically determine Q (t) in section IIIII 



II. CALCULATING THE PROBABILITY 
DENSITY FUNCTION 

The statistical properties of an RTN process have 
been studied extensively in the context of reliability 
theory, alternating renewal processes and queueing the- 
ory [ii, m, m, [li, [13, [li- In our case, we are spccifl- 
cally interested in the probability of the RTN spending a 
certain fraction of the observation period in a particular 
state. The PDF for an RTN signal fluctuating between 
the states and -1-1 is given by [23 as 



p (r, T) = Ae 



-AT 



T-T 



h 



(18) 



where Ii is the modifled Bessel function of the flrst kind. 
This PDF assumes the initial state is -1-1, and that at least 
a single fluctuation occurs. Here, r is used to describe 
the time spent in the state and T is the duration of 
the process we are considering. The parameter A charac- 
terises the fluctuator rate as before. Properly normalised 
the PDF is. 



fc=i 



-h 



2\^T {T~t) 



T 



T sinh^ {^) 



(19) 
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We could equally describe a process which begins in the 
state 0, with T — t describing the time spent in this state. 
Assuming at least a single fluctuation occurs, the full 
PDF is obtained by averaging over both possible starting 
states 



p' {t,T) = -\p{t,T)+p{T-t,T)]. 



(20) 



We may re-express this PDF in terms of the mean fluc- 
tuator state ^, where ^ S [— Ij !]• Taking care to preserve 
the normalisation, the PDF for an RTN process of dura- 
tion T assuming at least one fluctuation occurs is 



2^ 



AT 



T 

2«- 



yr^sinh^ (^) 



(21) 



(22) 



The case where no fluctuations occur must be treated 
separately. In this case we expect ^ to be either of ±1. 
The properly normalised PDF for this case can be de- 
scribed using two delta functions, 



^0 (e,T) = i[<5(e 



1) 



-^(e + 1)] 



(23) 



The full, general PDF is constructed by appropriately 
weighting fio (^,T) and il>o (C,?^), with the fluctuation 
probability given by the Poisson distribution 



Pk (AT) 



-AT 



(AT)* 



fc! 



(24) 



where k denotes the number of fluctuations, such that 

n it T) - po (AT) f7o (e, T) + p>o (AT) r!>o (e, T) , (25) 

and p:^o{XT) = 1 — po (AT). After simplification, the 
resulting PDF is 



-XT 



[<5(e-i) + ^(e + i)] 



AT 



1 



(26) 



In what follows we examine the three limiting cases 
of the PDF and use these to construct an approximate 
PDF. The approximate PDF provides greater physical 
insight when working within these limits. 

Examining the two limiting cases of the PDF 
^>o (CjT), the fast and slow fluctuator limits and com- 
bining them with fio (^, T), leads to a simplified expres- 
sion which approximates f2(^,T). We begin by consid- 
ering the slow fluctuator limit A ^ for the distribution 
describing at least one fluctuation, ri>o(^,T). This is 
the regime where no more than one fluctuation occurs. 
In this limit 



laix) 



1 



r(a + 1) V2 



(27) 



and 

sinh {x) =x + (x^) 
This reduces the PDF to 



^>o (^,T) 



1 

2' 

f^i(e,r). 



(28) 

(29) 
(30) 



This uniform distribution implies that a fluctuation is 
just as likely to occur at any time during the system 
evolution. 

The limit A — > oo represents a fast fluctuator. In this 
regime 



and 

sinh (x) 
which reduces the PDF to 



(31) 
(32) 

(33) 
(34) 



In this limit ^ will be small, therefore making the sub- 
stitution ^ ~ 1/ \J AT, we find the PDF to be Gaussian 
about the origin. 




f^>o (C,T) 



1 



n>i(^,T) 



(35) 
(36) 



which we expect intuitively. We note that this approach 
is similar to that used by Happer and Tam when consid- 
ering the Gaussian limit of rapid spin exchange in alkali 
vapors [29j . 

Weighting these two limiting cases and the PDF de- 
scribing no fluctuations using the Poisson distribution as 
before, allows us to construct an approximate PDF 

^ (C, T) ~ Po (AT) r!o (e, T) + pi (AT) (C, T) 

+ p>i(AT)f^>i(e,T), (37) 

where p>i (AT) = I - po (AT) - pi (AT). This ap- 
proximate f2 (^, T) provides nice analytic solutions for 
Q (<) in each of the three interesting fluctuator regimes. 
While this is only an approximation to the exact solu- 
tion (Eq. I26p . it can provide more physical insight, as 
will become apparent later. 
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III. USING THE PDF TO DETERMINE Q (t) 

The superoperator Q (t) can be derived analytically 
via Eq. [T7] using the PDFs determined in the previous 
section. Of particular interest is the non-unitary part of 
the superoperator. 

The non-unitary superoperator can be found by ex- 
panding the superoperator into a unitary and non- 
unitary part, such that Q (t) = Q^^^ (<) Q("^) (i). The 
evolution in the absence of a fluctuator is contained 
within the unitary part, 



Q(") (t) = e 



(38) 



whilst the effect of the charge fluctuator is contained 
within the non-unitary part Q^"") (i). Note that these 
two parts can be factored out due to the commutation 
relation (Eq. I12[) . We now determine the non- unitary 
parts of the superoperator for Q (^, T) and its various 
approximations. 

Beginning with the case where no fluctuations occur 
and the PDF is given by fig (C, T), as in Eq. [23l wc find 



(t) = cos (acriit) , 



(39) 



where cth is the Heisenberg superoperator introduced 
earlier. When there is at least one fluctuation (see Eq.| 
the resulting form of the superoperator is 



qIt^ it) 



{acT^tf - {\TY 



cos (acTfji) 



2 sinlr 



(40) 

Examining the limiting cases of the general PDF we find 
that for the slow fluctuator (see Eq. [50]) 



Q^' it) 



sin (acTfji) 
acrnt ' 



and in the fast fluctuator limit (see Eq. [36)) 



Q'>r (t) 



-(QM'J'Ht) /2 



(41) 



(42) 



It should be noted that this final superoperator corre- 
sponds exactly to that which would be obtained using 
the Lindbladian formalism [sO, [Hj if a decoherence oper- 
ator of the form L = "^""i ' '^2 was included. From this 
one can deduce that the fast fluctuator limit is equiva- 
lent to purely Markovian decoherence due to interaction 
with the environment via an exchange like two-qubit de- 
coherence channel. This is in contrast to conventional 
dephasing which is modelled using two independent az 
channels, one for each qubit. This distinction is particu- 
larly important as it implies that exchange fluctuations 
due to environmental charge fluctuations introduce cor- 
related errors which can have important implications for 
fault-tolerant QEC [H. 

Using the previous results, we can determine Q(™' (t) 
for the full weighted PDF's in both the approximate and 



exact cases. The exact PDF given in Eq. [26] yields 

2 



Q("^) (t) = e-^^ cos (acrnO 



(a^Ht)' - (AT)' 



- cos (acTHt) } , (43) 



whilst for the approximate PDF given in Eq. [37] we find 



Q("^' {t) 



-^^cos(a^Ht) + ATe-^^^^"("""^) 



l-e 



-XT 



- \Te-^' ) 



-{atiCrntf/2 



(44) 



In general it is difficult to graphically compare these an- 
alytic forms of the superoperator to the numerical result. 
However, it is possible in this CcLS6, 3jS the superopera- 
tor Q(™) (i) is a sparse matrix with the same underlying 
structure of the cth superoperator which defines it. It fol- 
lows from the definition of cth (see Eq. [6]) that the only 
non-zero matrix elements of the Heisenberg superoper- 
ator are ±2. Consequently, a comparison of the result- 
ing non-zero matrix clement of Q(""' (t), denoted Qnu, 
proves effective in determining the agreement between 
the analytic (exact and approximate) superoperators and 
exact numerical solution for the superoperator. 



AT = 100 



Z 




numerical o 
exact 
approximate 



2.5 



FIG. 1: Comparison of the analytical [exact (Eq. I43|l and 
approximate (Eq. I44p ] and exact numerical (simulated from 
Eq.[l6} solutions of Q'""^ (t). Plotted is the non-zero matrix 
element of Q^""-* (t), denoted QnU) as a function of time for a 
range of fluctuator rates, which span each of the three limit- 
ing regimes. The results show very good agreement between 
all three solutions, with the analytic approximation deviating 
only slightly when the fluctuator rate is on the time scale of 
T. 

The results, as a function of time for a range of fluctu- 
ator rates, shown in Fig. [H reveal very good agreement 
between the exact analytic and numerical results for all 
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rates A. The approximate solution also matches closely, 
particularly in the slow and fast fluctuator limits. Slight 
deviations from the exact solution can be seen when the 
fluctuations occur on the time scale of the process we are 
considering (AT « 1). In this regime the contribution 
from the uniform distribution fli (^, T) is at its maxi- 
mum and approximately on par with contributions from 
the other two distributions. The deviation from the ex- 
act results do not come as a surprise, as the approximate 
PDF is constructed from contributions due to 0, 1 or 
many fluctuations. Adding contributions from 2, 3 and 
more fluctuations would reduce this discrepancy. 

We now present the generalisation of the single fluc- 
tuator formalism to multiple fluctuators in the following 
section. 



ensuring that the spectral distribution function is prop- 
erly normalised 

J I S{a,\) dadX = 1. (49) 

The effects of a region of charge noise can now be mod- 
elled using either approach (Eq.|35]or Eg. BS)) . The choice 
will depend on exactly what information is known about 
the system. 

In the following section we outline how to implement 
the superoperators and provide a discussion of some im- 
portant limitations of this approach. 

V. USING THE SUPEROPERATORS 



IV. MULTIPLE FLUCTUATORS 

Extending this formalism to multiple fluctuators is 
straight forward and provides a method for the treat- 
ment of many physically realistic scenarios. The total 
ensemble averaged superoperator, A{t), for N fluctua- 
tors is just the product of all the individual ensemble 
averaged superoperators, Q (t), such that 



N 



A(i)=Q(") (onQi""^ w 



(45) 



i=i 



This result is useful for a finite number of fluctuators 
each with known strength and rate. However, in most 
instances only the spectral distribution in strength and 
rate will be known and therefore Eq. offers no fur- 
ther insight. By considering all possible unique super- 
operators Q("^' (tti, Ai, t) weighted by their probabihty 
of occurrence pi (where pi S [0, 1]), in a similar way to 
the method used to construct Q {t) in Eq. [iTl we may 
re-express Eq. US] as 



M 

A(t) = Q(") (i)n[Q'"'^ {a^,X^,t) 



i=l 



Npi 



(46) 



where in general there are M possible fluctuator types 
and N fluctuators. We would like to interpret the pi as a 
spectral distribution function in ai and Xi . By expressing 
A as a sum of logarithms 



M 



A (t) = Q(") (t) cxp <NY,P^ In [Q'""^ {a^,K, t) 

(47) 

and extending the definition of A to the continuum, we 
may replace the pi with a spectral distribution function 
S (a. A) such that 

A(<) = Q(-) (i) 
xcxpiN ff S (a. A) In [q("") (a. A, t)] da dx\ , (48) 



The superoperators in sections IIIII and IIVI were con- 
structed on the basis of the commutation relation for 
our Heisenberg Hamiltonian (see Eq. [T^ . The commuta- 
tion relation meant we could express the total evolution 
operator as a product of two evolution operators, each 
describing the total time spent in one of the fluctuator 
states (see Eq. [T3]). The superoperators themselves also 
commute as a result. Without further approximation the 
superoperators can be used to model individual processes 
satisfied by the Hamiltonian in Eq. [T^] (or a similar com- 
muting Hamiltonian, such as the Ising interaction). 

In more complex superoperator applications, for ex- 
ample when multiple applications of the superoperators 
are separated by non-commuting operations, it may be 
necessary to make a further approximation. When mod- 
elling these more complex processes (see Fig. [2] for an 
example where two superoperators are separated by a 
non-commuting gate operation, G) a problem arises with 
the formalism when considering the slow fluctuator limit. 
Specifically, the superoperator terms which should de- 




FIG. 2: Multiple instances of the superoperator Q (t) sepa- 
rated by a non-commuting gate operation. Attempting to use 
the superoperators to determine the effects of charge noise in 
a process like this can lead to the introduction of errors in the 
slow fluctuator regime. 

scribe no fluctuations at all, actually account for the pos- 
sibility of a fluctuation occurring between superoperator 
applications. We refer to these terms as cross terms. As 
the fluctuation rate increases the Poissonian weighting 
of these cross terms in the overall superoperator reduces, 
hence reducing the cross terms significance. We now con- 
sider a simple example which illustrates how these cross 
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terms manifest themselves, before showing how an ap- 
proximate solution can be constructed for the slow fluc- 
tuator regime by removing the cross terms. 

When the supcropcrators do not commute as in the 
example shown in Fig. [21 two or more applications of the 
supcropcrators leads to the introduction of unphysical 
cross terms in the slow fluctuator limit. This becomes 
apparent when we consider the action of the superoper- 
ator describing no fluctuations Qq""^ (t), which may be 
expanded in terms of the supcropcrators, (t), each 
describing one of the two fluctuator states in the no 
fluctuator limit 



Q(nu) 




W = I [QJ (t) + Qo (t)] ■ (50) 



The cross terms from the product of two (or more) of 
these supcropcrators which sandwich non-commuting op- 
erations results in the description of a single (or multi- 
ple) fluctuation(s). For example, consider the system in 
Fig. [21 The gate operation G does not commute with the 
supcropcrators 



Q 



total 



Q(t2)GQ(ii), 



(51) 



where G is the superoperator representation of the gate 
operation G. Expanding out each of the fluctuator su- 
pcropcrators using 



Q(t) = Q(") {t)[po (At)Q[,"") it) 

+ P>o(At)Ql:™^ (t) 
and Eq. 1501 with some rearranging we find 



(52) 



Q 



total 



Pa [A (ti 



*2)]r^(«) 



Q^"' (^2) [Q([ {t2) GQ+ (ti) 

Q+ {t2) GQo (ti) + Qq- ih) GQ+ ih) 

+ Qo ih) GQo ih) ] Q(«) (ti) + . . . , (53) 



where we have only shown the terms which should de- 
scribe no fluctuations. This entire expression should rep- 
resent the total superoperator describing no fluctuations. 
However, careful inspection shows the presence of two 
cross terms, which actually imply the occurrence of a fluc- 
tuation during the non-commuting gate operation. Cross 
terms of this form are actually a manifestation of this 
superoperator formalism and should be removed without 
also removing any unitary evolution. 

It should be emphasised that this problem only occurs 
in the slow fluctuator limit, where there is a signiflcant 
probability of there being no fluctuations during a two- 
qubit operation. As the fluctuation rate increases, the 
probability of a fluctuation occurring during the SINGLE 
qubit gate increases, which means that each application 
of the two-qubit gate becomes statistically independent. 
In this limit, the formalism as presented so far is exact 
and does not require any attention to cross terms. 



We now wish to remove cross terms describing pro- 
cesses which should not occur, such as those in Eq. [531 
When two or more applications of the Qq""^ (<) superop- 
erator occur in succession, the introduction of these un- 
physical terms also occurs. It is possible to construct an 
approximate solution by carefully removing these cross 
terms. 

Cross terms in the large A limit do not pose a prob- 
lem as the Q>'o ' (t) superoperator provides the domi- 
nant contribution to Qtotal in this limit. It will there- 
fore be most important to remove the cross terms due 

solely to the Qq""^ (t) superoperator - the zeroth order 
cross terms. In general, when many applications of the 
supcropcrators are required, there will be higher order 
cross terms. For example, the first order cross terms 
would contain the Qq'^"'' (t) superoperator and a single 
instance of the Q^g ^ (t) superoperator. Successive im- 
provements to the approximate solution arc achieved by 
removing these higher order cross terms. As more orders 
of these terms are removed the approximation improves 
for larger A, with the actual region of improvement de- 
pendent on the Poissonian weighting of the cross terms 
being removed. 

As an example, the zeroth order cross terms, denoted 
Xq, are removed by firstly re- weighting Qtotal using the 
Poisson distribution over the total process duration Ttotai 
(i.e., including non-commuting operations). 



Q 



total 



PO (ATtotal) Qtotal + P>0 (ATtotal) Q 



total- 



(54) 



The cross terms can now be removed 



Qtotal ~ PO (ATtotal) (Qtotal — Xq) 

+ P>0 (ATtotal) Qtotal, (55) 

taking care to not also remove any unitary evolution. The 
cross terms may also be removed from the second term, 
however the improvement from doing this is minimal due 
to the small contribution from the zeroth order terms at 
large A. 



VI. CONCLUSION 

The exchange interaction is of fundamental importance 
for controllable quantum mechanics in solid-state sys- 
tems. Its application to mediate spin flips or entangle- 
ment has particular importance in spintronics and QIP, 
hence the stability of the exchange interaction is crucial 
for precise time varying control. In solid-state spin sys- 
tems this stability can be affected by the local charge en- 
vironment, in particular charge fluctuators, due to the ex- 
change couplings dependence on the electron wave func- 
tion overlap. 

We have developed a model to describe the effect of 
charge fluctuators on the exchange interaction as a func- 
tion of time, using supcropcrators dependent on the noise 
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amplitude and rate. These superoperators can be in- 
cluded in time-dependent calculations of the state of the 
spin system to model the effect of the charge noise. Fur- 
thermore this analysis holds for other spin couplings, like 
the Ising interaction, where a commutation relation anal- 
ogous to Eg. [T2I exists. 

In the fast fluctuator limit we demonstrated how inter- 
action with the environment via an exchange like deco- 
hcrence channel leads to purely Markovian decoherence, 
although the decoherence operator leads to correlated 
noise across the two spins. 

The generalisation to multiple fluctuators means that 
the effect of charge fluctuators distributed according to 
a spectral distribution function can also be modelled. In 
the simpler case where only a small number of well de- 
fined fluctuators exist, the total superoperator is just the 
product of the individual fluctuator superoperators. 

As our model is completely analytic, the effects of ex- 



change fluctuations can in most instances be included 
trivially in more sophisticated analyses, without the need 
to explicitly sum over noise histories. This is important 
for analysing the operation of spintronic devices as well 
as QEC and fault-tolerance for QIP. 
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